Radiative capture of polarized neutrons by polarized protons 
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A model-independent irreducible tensor approach to pifi, 'y)d is presented and an explicit form for 
the spin-structure of the matrix M for the reaction is obtained in terms of the Pauli spin-matrices 
cr(n) and cr{p). Expressing the multipole amplitudes in terms of the triplet — * triplet and singlet 
triplet transitions, we point out how the initial singlet and triplet contributions to the differential 
cross section can be determined empirically. 
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The study of radiative capture of neutrons by protons 
is of interest not only to nuclear physics as a testing 
ground for theories of NN interaction, but also to astro- 
physics, where the fusion reaction is part of the proton- 
proton chain responsible for the generation of solar en- 
ergy and production of elements in the early universe 
101 . The cross section for the process was directly mea- 
sured 1^ for the first time at neutron energy of 550 keV, 
although measurements j^] with thermal neutrons have 
been carried out earlier. The 10% discrepancy noted be- 
tween such measurements and theory was accounted for 
with surprising accuracy by the inclusion |^ of meson ex- 
change currents (MEG) . Measurements at higher ener- 
gies were used to test effects of p and uj exchanges and rel- 
ativistic corrections to the impulse approximation |^ . Al- 
though a thorough review of the inverse reaction fi(7, n)p 
over a wide range of energies is found in Q , the energy re- 
gion just above threshold does not seem to have received 
much attention. Photodisintcgration experiments at low 
energies focused attention |^ on the relative Ml and El 
contributions. Experiments employing polarized photons 
have been reported between 5 to 10 MeV and more 
recently at 3.58 MeV. The angular distribution as 
well as polarization of the neutron in d{'-f,p)n were mea- 
sured pl| . Though the measured angular distribution of 
the neutron was found 1 12 to be in good agreement with 
theory at 2.75 MeV, the angular distribution and neutron 
polarization in ||l^ differ from theory which includes the 
MEG contributions. Employing polarized neutrons at 6 
MeV and 13.4 MeV, the analyzing powers in p{n, j)d 
was measured which differed from theory. Recent cross 
section calculations ^3, 14 1, which agree with each other 
within 5% deviation, are found to differ from the 1967 
estimates of Fowler et al ||l| . The cross section was ob- 
tained in by fitting the existing data with a poly- 
nomial expansion, while the calculation in 1 14 includes 
MEG's, isobar currents and pair currents. The theory is 
in good agreement with cross section data for neutrons 



TABLE I: Transitions from initial states of np system, from 
threshold onwards, to the final state with spin-parity 1"*" 
and isospin / = 0, together with the corresponding A7, As 
and multipole characteristics of the radiation emitted in the 
fusion reaction. 
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AI 


As 


Allowed multipolarities 
of emitted radiation 


'So 1 = 1 
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Ml 


^Si 1 = 








Ml E2 


'Pi 1 = 
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El M2 


'Po 1=1 
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El 


^Pi / = 1 
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El M2 


^P2 1=1 
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El A/2 E3 
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and so on. 



with energy above 14 MeV, but deviates by about 15% 
from older d{^,p)n measurements between 2.5 MeV and 
2.75 MeV pa, which correspond respectively to neutron 
energies 550 keV and 1080 keV in p{n,^)d. Apart from 
the fact that uncertainties in p{n,j)d cross sections at 
energies below about 600 keV lead to dominant uncer- 
tainties in the determination of the relative abundances 
of elements in the early universe, it is interesting to study 
the fusion reaction for its own sake, especially since a be- 
ginning has also been made H] to study p{n, ^)d precisely 
in the region of energies of interest to astrophysics. This 
recent study Q has shown experimentally that the Ml 
transition from S-wave capture, which dominates at lower 
energies and the El transition from P-wave capture at 
higher energies are comparable at 550 keV. 

It is interesting to observe that the neutron capture 
in p(n, 7)0? from the '^S'l state is from the isospin / = 
state, which is the same for the deuteron, while that from 
the ^^o, I — 1 state is characterized by a A/ = 1 transi- 
tion. Likewise, capture from ^Pi, I — state leads to a 
AI = transition, whereas those from ■^Pj=o,i,2, I — 1 
states are characterized by AI = 1 transitions (see TA- 
BLE m. The cross section at very low energies is dom- 
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inated by the A/ = 1 amphtude for the Ml transition 
from the ^.^o state, since the NN scattering length in 
the ^5*0 channel is large and the Ml amplitude is propor- 
tional to the isovector magnetic moment of the nucleon, 
which is more than five times larger than the isoscalar 
magnetic moment. It was, however, noted quite early 
| p6[ that the spin-dependent effects in the fusion reac- 
tion are sensitive to the small AI = amplitudes Ml 
and E2. Hence a "polarized-target-beam test" was pro- 
posed for the '^Si Si radiative transitions in ther- 
mal np capture. The isoscalar Ml and E2 amplitudes 
have also been studied theoretically recently using 
effective field theory, since the circular polarization Q 
of photons emitted in the capture of polarized neutrons 
by unpolarized protons and the angular distribution ^] 
in the capture of polarized neutrons by polarized pro- 
tons are sensitive to the presence of these amplitudes. 
Parameter-free predictions, employing the Weinberg | |2C| ] 
scheme of power counting, have also been made for 
the spin-dependent observables in p(n, 7)d at threshold 
energies. However, the recently measured value of 
rj = (1.0± 2.5) X 10~^ for 7-anisotropy at 50.5% polariza- 
tion of neutrons and protons in p(n, 'y)d cannot test the 
effective field theory predictions at 10"''. Although this 
measurement provides the first experimental value for a 
spin-dependent observable in p{n,j)d, spin-observables 
in elastic NN scattering have been measured for more 
than two decades More recently, attention was fo- 
cused on such observables in NN scattering, with a veiw 
to determine the exact strength of the tensor interaction 
and the case of 7r° production in pp collisions for 
which the threshold itself is high. Even at higher en- 
ergies, it is clear that the singlet and triplet radiative 
captures from any arbitrary initial partial wave, £ lead 
to pin, "f)d characterized respectively by AI = and 
AI = 1 or AI = 1 and AI = depending on whether the 
initial parity (—1)^ is odd or even. Since the NN inter- 
action, as elicited from elastic scattering data, conserves 
channel spin, s (apart from total angular momentum j) 
as well as isospin /, it would be of interest to study ex- 
perimentally the relative strengths of the initial singlet 
and triplet contributions to the p{n,j)d reaction at any 
given energy, since the fusion reaction leads to transitions 
in the two-nucleon system, which change spin s as well 
as isospin /. It is also worth noting that np fusion near 
threshold bears some resemblance to the case of pd fu- 
sion, where the validity of the so-called "no-quartet the- 
orem" Q was questioned and led to several incisive the- 
oretical studies [g6[ in later years. More recently, it was 
pointed out p7| that a non-zero tensor analyzing power 
in p {d ,^)^He by itself provides a clear signature to 
the contributions from the quartet amplitudes. More- 
over, it was shown that it is possible to determine 
empirically the individual doublet and quartet differen- 
tial cross sections at any given energy through appropri- 
ate measurements of the relevant spin-observables in Nd 
fusion. In contrast to the difficult Nd fusion experiments 



suggested in it should be technically more feasible 
to carry out the p{n, "f)^II experiments Q. 

The purpose of this Rapid Communication is to present 
a model-independent theoretical approach to the np fu- 
sion process based on the irreducible tensor formalism 
and identify observables in p{n, d)j experiments that fa- 
cilitate determination of the singlet and triplet cross sec- 
tions empirically for the fusion reaction at any energy. 

Let p{p, 0, 0) and k{k, 9, 0) denote respectively the neu- 
tron and photon momenta in a right-handed cm. frame 
p9[ and let (md; fc/i|r|p; m„, rrip) denote the elements 
of the on-energy-shell T-matrix for the fusion reaction 
p{n^j)d where mn,mp,md denote respectively the spin- 
projections of the neutron, proton and deuteron and 
II — ±1 denote the left- and right-circular states of pho- 
ton polarization, as defined by Rose pof . The unpolar- 
ized differential cross section for the reaction is then given 
by 



do-Q 



E„EpEd 1 



\2tiE ) p 4 



E E 



(1) 



mn,mp)\ , 



where Ep,En and Ed denote the cm. energies of the 
proton, neutron and "H respectively, when the reaction 
takes place at cm. energy E. Expressing the initial np 
system in terms of isosinglet and isotriplct states and 
making use of the standard multipole expansion | ]30| for 
the photon in the final state and the usual partial wave 
expansion for relative motion in the initial state with 
channel-spin s = 0, 1, the T-matrix may be expressed as 

^ = E E E {S\l,s)-T\,j„s)), (2) 

s=Q,l A=|l-s| M = -l,l 

where S^{l,s) denote |^ irreducible tensor operators 
of rank A in hadron spin-space. The irreducible tensor 
amplitudes T^in, s) in Eq. (|^) are given explicitly by 



x(-l) 



00 00 

(2-)(-i)^-^E E E 
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where the energy dependence is carried entirely by the 
partial wave magnetic and electric 2^-multipole ampli- 



tudes T 



3 (mag) 
L-ls 



(E) and Tij'f°'^'^ (E) respectively, while the 



angular dependence is contained entirely in d^p^{9). We 

use the shorthand [j] = ^ {2j + 1) and the rest of the no- 
tations follow pO|. The multipole amplitudes are given 
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in terms of the reduced matrix elements by 

Ti;r\E)^^-{{-l)^+' + {~lY} {\\T\\) 



1 



+ (-!/} (||T||), 
where ( ||r|| ) denotes 

(11^11)= (-i/w-^ 

7=0,1 

X C(ii/;-iiO)((lL)j;/d = 0||r||(^s)j;/ 



(4) 



(5) 



Expressing the irreducible tensor operators S^(\,s) of 
rank A in terms of the Pauli spin-matrices cr{n) and (t{p) 
and unit 2x2 matrices crQ(n), ctq{p) for neutron, proton 
respectively following ||3l|] , we may rewrite Eq. (0) in the 
elegant form 



^.iTr(MMt 
dil 4 V 



(6) 



where Tr(= X^m^ '^n) denotes the Trace or Spur and the 
matrix M for has the form 



1 (A1+A2) 

Ai,A2=0 A=|Ai-A2| 



(7) 



in terms of irreducible tensor amplitudes Af^(Ai, A2; /i) 
of rank A given by 



M,^(Ai,A2;/i) 




where { } denote Wigner-9j' symbols js^]- We may ex- 
plicit Eq. ^ as 



M =A + B (cr(n) • a{p)) + (cr(n) + ct{p)) ■ C 
+ (cr(n) - cr(p)) • D + {cr{n) X cr{p)) ■ E 

+ ({a{n)®CT{p)f -F^), 

where the coefficients are related to (|^) through 

1 



(9) 



Mo"(0,0;m) ; B 



1 



V3 



Mo(i,1;m) 



(10) 



Cl^^ [m1(1,0;m) + M1(0,1;m)] ; 
^i = ^K'(l,0;M)-Mi(0,l;/i)]; 



andC3, -D^, Sj'; denote respectively the spherical compo- 
nents of C, £>, E. Comparing Eq. (^) with Af for elastic 
NN scattering shows clearly that the fourth and fifth 
terms containing {cr{'n) — cr{p)) and {cr(n) x (t(p)) are 



the ones which induce transitions from the initial spin- 
singlet state to the final spin-triplet state of the deuteron. 
To estimate the singlet and triplet cross sections empir- 
ically we observe that the differential cross section for 
p(n, 7)0? is given by 

^ = Tr(MpMt), (11) 

where the density matrix 

p=i[l+(^(n).P(n))] [l+{a{p)-P{p))] (12) 

describes the initial spin state if P{n) and P{p) de- 
note respectively the neutron and proton polarizations. 
Rewriting Eq. (01]) as 



(13) 



E Po.{n)Pp{p)B^p, 

a^l3—0,x,y,z 

where Po{n) = 1, Po(p) = 1 and 

B^p = Tr (Ma„(n)a0(p)M^) , (14) 

we readily see that the unpolarized differential cross sec- 
tion d) is given by (l/4)Boo- Noting pi further that 



7r(l, 1) = 1 [1 + a,{n) + a,{p) + a^{n)(j,{p)] 



(15) 



7r(l,0) = ;^ [1 + (^x{n)(Tx{p) + ay{n)ay{p) ~ az{n)az{p)] 

(16) 

7r(l, -1) - J [1 - (^z{n) - <Jz{p) + <7,{n)a,{p)] (17) 



7r(0, 0) = - [1 - a.j;{n)a^{p) - (jy{n)(Ty{p) ~ a^{n)a^{p)] 



(18) 
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are the projection operators \sm){sm\, s = 0, 1; m 
+s, . . . , — s, we readily identify 
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(22) 



as the triplet and singlet contributions (d(Ts.„i)/(df2) 
which add up to (dcro)/(df7) given by Eq. (0). Eq. (0) 
can also be expressed in the form 

do_ 

dn ^ ^ ^ 



' dn 



J2 ({P''-{n)®P^^{p)f -B^K.k,,)) (23) 

k 

1 + P(n) • A{n) + P(p) • A{p) 
+ ^2 ({P{n)^Pip)f-A'^] , (24) 



A;=0 



where the irreducible bilinear amplitudes in Eq. 



X E (-l)^WM[A][A']W^(s'fclA;sA') 

s.s' , A, A' 
1 1 

fc 



2 2 



X 
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(T\/i,s)®rt''(^,s') 



(25) 



are related to the spherical components of the analyzing 
powers in Eq. (p3) through 



^Aj(n) =5^(1,0) 



(26) 
(27) 
(28) 



with i?Q(0, 0) = (d(To)/(dri), the unpolarized differential 
cross section. From Eq. (^), we observe that the T^ifJ-, s) 
satisfy the symmetry property 



r,^(A*,s) = (-i)^-^r^,(-M,s) 



(29) 



This in turn implies through Eq. ( p5| ) that the bilinear 
amplitudes satisfy the property 



B„ {km kp) — ( 1) B_ {kn, kp) 



(30) 



so that 

^0 {kn, kp") — 
and from Eqs. (|^) and ([27|), 

A,{n) = A,{p) = 0. 



(31) 



(32) 



Using these results, we can now write the channel-spin 
differential cross sections (di7s_m)/(dri) in terms of the 
analyzing powers as 



dcri,i 


d(Ti,-i 


dn 


dn 


dcri,o 


_ 1 ddo " 


dn 


~ 4 dn . 


dcro,0 


_ 1 do-Q " 


dn 


~ A dn . 



_ 1 dffo r 
^ A dn 

l + \/3A° 



1 



J_/lO 



3^0 



2 42 



(33) 
(34) 
(35) 



The analyzing powers, moreover, are measurable readily 
in an experimental setup such as in ||l9|| , where P{n) = 
P{p) = PP^ with P ^ 50.5% and P could be longitu- 
dinal as well as transverse. With P chosen parallel and 
antiparallel to the beam, it is clear that 



da 



dn 



do--z 

dn 



dn 



1 



V3 



{V2AI-AI) 



(36) 



which yields (v^^g - Ajj). 

Likewise, if P is chosen parallel and antiparallel to a 
direction (say x of y) perpendicular to the beam, we have 



dcr+:i 

dn 



dn 



A(J-x 

dn 



dn 



p2 



(37) 



dq--y ^ 2 — 



dn 



dn 



1 - 



p 



(38) 



where we have made use of the fact that A2 = At 



, from 

Eqs. (g^) and iM). By adding Eqs. (|7|) and one 
readily obtains (^{] + ^^o)- Since (v2A^ - ^0) is al- 
ready known from Eq. (p6[), one can determine Aq as well 
as Hence all the (dqs,m)/(dr2) for m — —s, . . . ,s and 
s = 0, 1 are determinable from experiment empirically at 
any given energy. We may perhaps add that these mea- 
surements need not have to be carried out at accuracies of 
10~^ to determine {das^m)/ (dn). Moreover, the discus- 
sion presented here aims to extend the program of Breit 
and Rustgi | |T^ to actual individual determinations of the 
three triplet and one singlet cross sections for the impor- 
tant 10] fusion process at the differential level itself rather 
that estimate the relative importance of the triplet am- 
plitude vis-a-vis the singlet amplitude. It would, there- 
fore, be desirable to extend the recent experiment jl^ on 
p{n, 'y)d to measure the observables (pq) to (p8|) in order 
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to determine the differential cross sections {das,m)/i<i^) 
individually and hence study the fusion reaction more 
incisively at any given energy. 

It is perhaps not out of place to mention here that the 
highly successful NN potential models like Nijmegenl 



and Nijmegen2 potentials 
and the Bonn potential [37 



35| , the Argonne potential |36f| 
which reproduce the elas- 
tic NN scattering data with a remarkable /datum w 
1, have been subjected to an interesting study of the 
deuteron properties by Polls et al. |38|, who say, It is well 
known that the off-shell behavior of NN potentials can- 
not be pinned down by on-shell data... Moreover, within a 
given model, pinning down the on-shell point limits the 
range of variation for the off-shell behavior. The four 
modern high-precision models considered in this study pin 
down the on-shell T -matrix as much as by all means pos- 
sible since they fit the NN scattering data with the per- 
fect x^/ datum « 1. Furthermore, it may be reasonable 
to believe that the four models cover about the range of 
diversity that there is to realistic physical models for the 
NN interaction. Based upon these premises, one may 
then conclude that the off-shell uncertainties revealed in 



this study are what we have to deal with, at the current 
status of theoretical nuclear physics. At this time, we do 
not know of any other objectively verifiable aspects that 
could further reduce the off-shell uncertainties. Tighter 
constraints for the off-shell behavior of NN may emerge 
in the future... In this context, it is perhaps pertinent to 
point out that np ^7 involves necessarily off-shell NN 
interactions. Therefore it would also be of interest to ex- 
amine the fusion reaction with a view to throw light on 
the off-shell uncertainties. Further work is in progress. 
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